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We analysed the onset temperature Tp for the pairing in cuprate superconductors at small dop- 
ing, when tendency towards antiferromagnetism is strong. We considered the model of Moon and 
Sachdev (MS), which assumes that electron and hole pockets survive in a paramagnetic phase. 
Within this model, the pairing between fermions is mediated by a gauge boson, whose propagator 
remains massless in a paramagnet. We related the MS model to a generic 7— model of quantum- 
critical pairing with the pairing kernel A(f2„) ex 1/0^. We showed that, over some range of param- 
eters, the MS model is equivalent to 7 = 1/3-model (A(r2) tx fi"^'^). We found, however, that the 
parameter range where this analogy works is bounded on both ends. At larger deviations from a 
magnetic phase, the MS model becomes equivalent to 7 model with varying 7 > 1/3, whose value 
depends on the distance to a magnetic transition and approaches 7 = 1 deep in a paramagnetic 
phase. Very near the transition, the MS model becomes equivalent to 7 model with varying 7 < 1/3. 
Right at the magnetic QCP, the MS model is equivalent to the model with X{Qn) oc logJl„, which 
is the model for color superconductivity. Using this analogy, we verified the formula for Tc derived 
for color superconductivity. 

PACS numbers: 



I. INTRODUCTION 

A quest for understanding of the phase diagram of 
cuprates and other strongly correlated electron systems 
generated strong interest to the pairing problem near a 
quantum critical point (QCP). In distinction to a conven- 
tional BCS/Eliashberg theory, in which pairing is medi- 
ated by an interaction which can be approximated by 
a constant at small frequencies, quantum-critical pair- 
ing is mediated by gapless, dynamical collective modes^. 
In most cases, only a particular charge or spin mode 
becomes gapless at criticality, and the same interaction 
that gives rise to pairing also accounts for the fermionic 
self-energy. In systems with d < 3, d'S{k,uj)/duj di- 
verges at oj — T' 0, i.e., interaction with a gapless collective 
modes destroys fermionic coherence either at particular 
hot spots ^, if the energing order parameter has a finite 
momentum q, or along the whole Fermi surface (FS) if 
the ordering is with q = (Refi^). Then, in another 
distinction to a conventional BCS/Eliashberg theory, the 
pairing involves incoherent fermions. 

The self-energy at a QCP generally behaves as uj^lo^'^, 
where 7 > 0, while the effective dynamic pairing inter- 
action (the analog of phonon a'^{n)F{fl)) scales X{il) oc 
{H.o/n)''. The kernel K{uj,uj') of the equation for the 
pairing vertex $(0;) — ttTJ^uj' K{uj,uj')^{uj') then scales 
as 



K{uj,uj') 



X{lu- 



cx 



-^/hl./li-T' 



(1) 



|S(a.')| 

It has the same scaling dimension —1 as in ECS the- 
ory, which corresponds to 7 = 0. However, in distinc- 
tion to BCS theory, where the kernel is l/|i^'|, the ker- 
nel if (a; — w') for 7 > depends on both external and 



internal energy, and the onset temperature for the pair- 
ing Tp has to be obtained by solving integral equation 
for the frequency-dependent pairing vertex $(a;) (we la- 
bel this temperature as Tp to distinguish it from the ac- 
tual superconducting Tc which is generally lower because 
of phase and amplitude fluctuations of $(w)). The fre- 
quency dependence of the vertex by itself is not unique 
to a QCP and is present already in the Eliashberg theory 
of conventional superconductors^'^. However, as we said, 
in conventional cases, the frequency dependence is rele- 
vant at high frequencies, while the low- frequency sector 
can still be treated within BCS theory. The new element 
of quantum-critical pairing is that non-trivial frequency 
dependence of the pairing vertex extends down to w = 
leaving no space for the BCS regime. 

To shorten the notations, below we label the model of 
QC pairing mediated by A(ri) oc l/|r2|''' as the 7-model. 
Examples of 7— models, include 2D pairing by ferromag- 
netic or antiferromagnetic spin fluctuations (7 ~ 1/3 
and 7 w 1/2, respectivelj^"— 1^, and 3D pairing by gap- 
less charge or spin fluctuations (7 = 0-1- , implying that 
X{n) ~ log 17 and E(w) ~ wlogw, Ref^ii^). Another ex- 
ample of 7 = 0-|- behavior is color superconductivity of 
quarks mediated by gluon exchange^. The 7=1/3 model 
describes the pairing of composite fermions at v — 1/2 
Landau leveljS, A similar problem with momentum inte- 
grals instead of frequency integrals have been considered 
in Ref. To". 

The solutions of the pairing problem in these systems 
differ in details but have one key feature in common - Tp 
monotonically increases upon approaching a QCP from 
a disordered side. At a QCP, Tp — fio f{j), where 7(7) 
is a smooth function, which decreases when 7 increases. 

Recently, Moon and Sachdev (MS) considered^ an- 



other example of QC pairing - the pairing of fermions 
in the presence of spin-density-wave (SDW) background, 
which can be either in the form of long-range SDW or- 
der, or SDW precursors. MS argued that the phase 
with SDW precursors is an algebraic charge-liquid in 
which fermions still posssess the same pocket FS as in 
the SDW phased. They found that the model remains 
critical away from critical point because of the pres- 
ence of gapless gauge fluctuations, and Tp actually in- 
creases as the system moves away from a QCP into a 
charge-liquid phase. This trend is opposite compared to 
a magnetically-mediated pairing without SDW precur- 
sors. This result is quite important for the understand- 
ing of the pairing in the cuprates as it shows that the 
onset temperature of the pairing instability does form a 
dome centered around a doping at which the system de- 
velops SDW precursors. i^ii^i^i In the region where SDW 
precursors are not yet formed, there is a large Fermi sur- 
face, and Tp increases as doping decreases and the system 
comes closer to a magnetic instability. Once SDW pre- 
cursors are formed above Tp, the trend changes, and Tp 
passes through a maximum and begins decreasing with 
decreasing doping. 

The reason for the opposite behavior of Tp in the pres- 
ence of SDW precursors is the suppression of the d— wave 
pairing vertex. In the ordered SDW state, the interac- 
tion between fermions mediated by Golstone bosons is 
dressed up by SDW coherence factors and vanishes at 
the ordering momentum by Adler principle^ In the dis- 
ordered phase, there are no such requirement, but the 
vertex is still suppressed as long as the system has SDW 
precursors. As the system moves away from a QCP, SDW 
precursors weaken, the pairing interaction gradually in- 
creases towards its value without precursors, and Tp in- 
creases. 

The goal of this communication is to relate the pairing 
problem in the presence of SDW precursors with earlier 
studies of the pairing within the 7 model. We argue that 
the MS model belongs to the class of 7-models, however 
the relation between the two is rather non-trivial. First, 
the pairing in the MS model remains quantum-critical 
even away from the QCP, when spin excitations acquire a 
finite mass ttt. at T = 0, i.e., the MS model away from the 
QCP is equivalent to a 7-model right at criticality. Sec- 
ond, we show that the MS model at different deviations 
from the QCP is equivalent to 7— models with different 7, 
ranging from 7 = 0-|- right at the QCP to 7 = 1 deep into 
the magnetically-disordered phase (but still with SDW 
precursors). Third, the scale fto in A(ri) = (rio/lf^l)''' 
also becomes 7-depcndent and increases with increasing 
7. This is the key reason why Tp ex Qq increases into the 
disordered state. 

We also consider in more detail the pairing problem 
right at the QCP. We find that at this point, the pairing 
interaction at small frequencies is X{Q) — g log{flo / \^\) 
and the fermionic self-energy S(u;) = gujlog{flo/\uj\), are 
logarithmical and depend separately on the cutoff scale 
r^o and on the dimensionless coupling constant g. This 



separate dependence is unique to the logarithmical form 
of the pairing interaction, for any power-law form of 
A(il) = (r^o/r^)''', the coupling g is incorporated into ujq. 
At frequencies H. > flo, A(ri) decays faster than the loga- 
rithm, eventially as l/fi. As we said, the same logarith- 
mical form of the pairing vertex appears in the problem 
of color superconductivity. For color superconductivity. 
Son obtained at weak coupling Tp ex ujQe"'^^'^^^^ [Ref.^-]. 
The result is similar to the BCS formula, but the depen- 
dence on the coupling is ^/g rather than 1/g. Son's result 
was confirmed in subsequent analytical studies, but, to 
the best of our knowledge, it has not been verified by 
actualy solving the linearized gap equation for the case 
when the interaction is logarithmical at small frequencies 
and decays as a power of l/ft at larger frequencies. We 
solved this equation for our A(f7), obtained Tp, and found 
excellent agreement with e^^/^^^-' behavior. 

The paper is organized as follows. In Sec. |lT]we briefly 
review (i) Eliashberg-type theory for the pairig induced 
by frequency dependent pairing kernel K(fl), (ii) the 7- 
model with no SDW precursors, and (iii) the MS model 
for the pairing in the presence of SDW precursors. In 
Sec. InT] we show that the MS model reduces to the set 
of 7— models with varying 7, which gradually increases 
from 7 = 0-1- as one moves away from the QCP. We argue 
that 7 remains 1/3 over some range of deviations from 
the QCP, but this range is quite narrow. In Sec. |IV] we 
consider in more detail the MS model at the QCP and 
show that at this point it is equivalent to the 7 model 
at 7 = 0-|- which in turn is equivalent to the model for 
color superconductivity. We present the numerical solu- 
tion for the pairig vertex, and show that logTc/Ep scales 
as 1/^- Sec. |V]prsents our conclusions. 



II. REVIEW OF THE MODELS 
A. Eliashberg theory 

Throughout the paper we assume that the pairing 
problem can be treated within Eliashberg-type theory, 
i.e., assume that the pairing kernel and the fermionic self- 
energy have singular frequency dependence but no sub- 
stantial momentum dependence. The justifications for 
neglecting the momentum dependence are to some extent 
problem-specific, but generally are due to the fact that 
collective bosons are Landau overdamped what makes 
them slow modes compared to electrons. 

Within Eliashberg theory, the fermionic Green's func- 
tion in Nambu spinor notations is given by 

E(w„) = iE(iu„)fo + $(w„)fi (2) 

where G^^ — Gq^ + S, f are the Pauli matrices in the 
particle-hole space, e^ is the fermionic dispersion in the 
normal state, $(u;„) is a pairing vertex, and Ti{ujn) is a 



conventional (non-Nambu) self-energy. The Nambu self- 
energy T,{lo„) in turn is expressed via the full Green's 
function as 

= ii;(w„)fo + $(w„)n 

The full set of Eliashberg equations is obtained by 
matching the coefficients of the Pauli matrices term by 
term: 



boson, and A(a;„ 



i;(a;„) = ttT^ A(a;„ -a;„) 
$(w„) = 7rr^A(w,„-i 



Urn + Y,{uJn) 



y/{uj^ + i:{uJmW + <^^{uJm) 



$(w™) 



(3) 

This set of two non-linear self-consistent integral equa- 
tions can be reduced to just one integral equation by 
introducing the inverse quasiparticle renormalization fac- 
tor Z(ujn) = 1 + S(cL'„)/a;„ and the pairing gap A(aj„) — 
$(w„)/Z(cj„). The equation for A(w„) then decouples 
from the equation on Z{ujn) and takes the form 



A(cj„)=7rT^ 



—====== A(a;„) - A(a;„) — 



(4) 
The inverse quasiparticle renormalization factor Z{ujn) is 
then obtained by substituting the result for A(a;m) into 



l;„Z(w„) = ^„ + ttT y^ A (cjm - iLj„ ) 



v/^^-f A2(^„) 



(5) 



To obtain Tp we will need to set $, A — > and solve 
the linearized equation either for the pairing vertex or 
the pairing gap. Both are non-self-consistent equations 
in this limit because at $ can be safely dropped from the 
expression for S(cl'„) which becomes: 



i;(w„) == ■nT'^\{uj„i~LOn)sign{uJrn), Z{ujn) = 1- 



(6) 
It is more straighforward to analyze the linearized equa- 
tion for the pairing vertex $(a;„): 

$(a;„)=7rT^ifK,cj„)$(a;„), (7) 

where the kernel is 

\uJrn\Z{uJm) 

B. the 7- model 

We use the term 7-model as abbreviation for a criti- 
cal model in which the pairing is mediated by a gapless 



has a power-law form 



a(o)h'^ 



(9) 



We assume in this work that 7 < 1. The analysis of the 
7— models with 7 > 1 is more involved (and more non- 
trivial) , but we will only need 7 < 1 for comparisons with 
the MS theory. 

To put this into perspective, we remind that in BCS 
theory A(r2„) is a constant up to a cutoff scale above 
which it is set to zero. In Eliashberg theory for non- 
critical superconductors A(ri„) has some frequency de- 
pendence but still reduces to a constant at small frequen- 
cies. In this situation, the pairing problem differs from 
BCS problem quantitatively but not qualitatively. In Eq. 
dl]), however, the pairing interaction preserve sigular fre- 
quency dependence down to J7 = 0. This makes the pair- 
ing problem in the 7— model qualitatively different from 
BCS problem. 

The divergence of A(f7„j) at zero bosonic Matsubara 
frequency 17 = is by itself not relevant for the pairing 
as can be straighforwardly seen from Eq. ^ for A(a;,„): 
the term with jn = n in the r.h.s. of this equation van- 
ishes. This vanishing is essentially the manifestation of 
the Anderson theorem for a dirty s— wave superconduc- 
tors because the scattering with zero frequency transfer 
is formally analogous to impurity scattering2i2^i^. The 
recipe then is to eliminate the term with zero bosonic 
Matsubara frequency from the frequency sums, what we 
will do. 

The fcrmionic self-energy for such \{n) has the form 



S(w„) 



no 



UJ„ 



Sh,n) 



(10) 



where 



S'(7,m) 



(m-Hl/2)i-T 



[C(7)-C(7,"^ + l)] (11) 



where (^(7) and C(7, 'tn + 1) are Rieman zeta function and 
generalized Rieman zeta function, respectively. We plot 
5(7, m) in Fig. [TJ At large m, S'(7, m) approaches 1/(1 — 
7), and the self-energy becomes E(w„) = aj^~'''(riQ/(l— 7) 
(this limiting form is reproduced if the summation over 
ujm is replaced by the integration). 

The rate with which S'(7,to) approaches 1/(1 — 7) by 
itself depends on 7 and slows down when 7 approaches 
one. For 1 — 7 << 1, S'(7,m) at large m behaves as 
5(7, m) w (1 — m~^^^'^^)/{l — 7) and becomes logm at 

7 = 1- ^ ^ 

Substituting the self-energy into the equation for 
^{ujjri) and rescaling the temperature by ilo we find, in- 
stead of d?!). 
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FIG. 1; Generalized Riemann-zeta function, S{'y,m). The 
dashed (red), dot-dashed (green) and dotted (blue) lines cor- 
respond to m = 100, 1000, 10000. The black thick line is the 
asymptotic function 1/(1 —7) valid at large m (see the text). 



It is a'priori not guaranteed that this equation has a so- 
lution at some Tp, but, if it does, Tp obviously scales 



as 



Tp = nofh) 



(13) 



Eq p2)l has been analyzed both analytically and numeri- 
cally^'^^^Sii^ii^, and the result is that (i) the solution ex- 
ists, and (ii) f{'j) monotonically decreases with increas- 
ing 7, and remains finite at 7 = 1, where /(I) = 0.254 
(RefiS). 



the MS model 
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FIG. 2: The momentum dependence of Db{k) at zero temper- 
ature (black line) very near the QCP (top panel, m — 0.01) 
and far away from the QCP (lower panel, m = 100) k is in 
units of 2kF and Db{k) is in units of m*v with ai = 0.5. The 
thin-dashed line (Red) and the thin-dotted line (Blue) show 
linear and quadratic behavior, respectively. At the smallest k, 
Dt{k) is quadratic for all rn > (see inserts), but for small m 
Dbik) raidly crosses over to the linear behavior. For large m, 
the quadratic dependence survives up to maximal k = 2kF- 
Note that the magnitude of Dbik) far away from the critical 
point is much smaller than near the QCP. 



The MS mode^ describes a phase transition between 
an SDW state and a magnetically disordered, alge- 
braic charge liquid which still preserves electron pockets 
near (0, tt) and symmetry related points in the Brillouin 
zone. The model is two-dimensional and involves non- 
relativistic fcrmions with small FS, relativistic bosons 
which describe antiferromagnetic magnons, and a gauge 
field which couples bosonic and fermionic fields. There 
are two non- relativistic fermions (two FSs), with the ac- 
tion 



Co 



A 



Or 



2m* 



M 



9±- 



(14) 



where m* is a fermionic mass, and the chemical potential 
/i is chosen to fix the Fermi momentum kp. 

The action of a relativistic boson (an antiferromagnetic 
magnon) is 



N 



a=l 



A=^ |9r^aP+«'|V 



z„|^ +m^|z„P 



where m is the mass of a magnon 



m = 2Tln 



4 + e — 



(15) 



(16) 



and the parameter a2 measures the distance to a QCP 
(we use the same notations as in^). For 02 < and T = 
0, the system is in the ordered SDW phase (m — 0), for 
a2 > 2 and T = 0, the system is in the disordered phase 
(m > 0). At T > 0, there is no long-range order, and 
TO > for all 02, although for 02 < it is exponentially 
small. 

The gauge field couples bosonic and fermionic fields 
via minimal coupling(9^ — > 9^ ± iA^). The two fermions 
{g±) have opposite gauge charges, and gauge- field in- 
duced coupling between fermions and magnons is attac- 
tive and give rise to pairing. Both longitudinal and trans- 
verse parts of the interaction contribute to the pairing, 
but non-trivial physics comes from the interaction me- 
diated by a transverse boson. For simplicity, we only 
consider transverse interaction. 

The form of the transverse propagator of the gauge 
field xtr{k,Q,n) is obtained by combining the Landau 
damping term coming from fermions with the static part 
coming from magnons. MS found, for the physical case 
of iV = 2 bosonic spinors. 



Xtr{k,Vtn) = 5.. 



Ixi-i nJo 



1 



2Db{k) + kF\nn\l{^k) 



(17) 



where 



Dbik) 



dx- 



1 



Stt Jo ^m? + v'^k'^x{l — a;) 
X coth 



^Jm? + v'^k'^x{l — x) 



2T 



(18) 



In Fig. [2] we plot Di,{k) at T = near and far away 
from a critical point. At the smallest fc, Di,{k) scales as 
fc^ for any tti > 0. For large m, this behavior holds up 
to k ^ kp, but at small m, the fc^ dependence only holds 
for the smallest k crosses over to Db{k) oc k at vk > m. 

The frequency-dependent pairing interaction A(ri„) is 
obtained in the usual way, by interating the gauge field 
propagator Xtr{k, ri„) in the direction along the FS: 



A(r!„) 



aiv 
2^ 






dk 



y/l - {k/2kFf 

2Db{k) + kF\^rM{^k) 



(19) 



where, following MS, we introduced another dimension- 
less parameter 



fcf 



ai 



(20) 



We will see that it plays the role of the dimensionless 
coupling. 

Because D},{k) interpolates between k and k^ depen- 
dences, the momentum integral in ([20l) diverges at fin — 
no matter whether the system is at the QCP or away 
from the QCP (i.e., whether m = or m > 0). This im- 
plies that the pairing in the MS model remains quantum- 
critical even when the system moves away from the actual 
critical point where SDW order emerges. Alternatively 
speaking, away from the QCP, magnetic excitations ac- 
quire a mass but the pairing gauge boson remains mass- 
less. 



III. THE MS MODEL AS THE 7 MODEL WITH 
VARYING EXPONENT 7 

As the pairing is generally a low-temperature/low- 
frequency phenomenon, it is instructive to analyze first 
the form of A(J7„ in the MS model, Eq. ((ID), at T = 
and in the limit of small r2„. Suppose that the sys- 
tem is in the paramagnetic phase, where m tends to a 
constant at T = 0. Using ([TC|. we obtain at T = 0, 
m = 477171* v'^a2 = 4:Trkpva2/ai. For small fin the mo- 
mentum integral in (|19l) is determined by k for which 
D\j{k)k/kF ~ \ft-n\- These typical k are small when 
i7„ are small. Substituting the small k expansion of 
Db{k) = v'^k'^/{Smn) (see Eq. ^) into ^ we obtain 



mn) 



n 





101 



1/3 



(21) 
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FIG. 3: The actual critical temperature Tp in the MS model 
with ai = 0.65 (the black line) and a\ ~ 0.75 (the dotted 
blue line) (from Ref.-^) and Tp''^^^^ from Eq. dgj) (the dashed 
red line) Although the trend is the same - both the actual Tp 



and r, 



7=1/3 



increase with increasing 02, the functional forms 



are quite different. The two functions just cross at 02 ^ 0.01. 
Besides, the actual Tp does depend on qi, particularly for 
small a2, while Tp^ is independent on ai (see the text). 



where 



n, 



7=1/3 




5127r2 
8IV3 



Efo^ 



(22) 



and Ef — kp/{2ni*). The corresponding self-energy is 



E(..™) = I u^y^ inr'^'y/\ 



(23) 



This X{ftm) is the same as in the 7-model with 7 = 1/3. 
The onset temperature for the pairing in the "1/3" model 

has been calculated in Refsi^ii^ and is Tp « 3.31S1q~ ' . 
Combing this result with Eq. (|22p . we obtain 

-2 (24) 



Tp^Tpl=^^^^ 119 Ef a^ 



Observe that T, 



7=1/3 



depends on ai only via Ef- 



Moon and Sachdev obtained Tp in the MS model nu- 
merically, using the full forms of Di,{k) and m{T). In 
Fig. [3] we compare their result with Eq. ([M|) . We see 
that while the trend is the same - Tp icreases with in- 
creasing a2, the functional forms of the actual Tp and 

Tp~ ' are very different, and Eq. ([M)) agrees with the 
numerical result for Tp only in a tiny range of 02 near 
a2 ~ 0.01. Besides, Tp/Ep obrained numerically does 
depend on ai, in distinction to Eq. ((M)) . 

We looked into this issue more carefully and found that 
the reason for the discrepancy is that Eq. (l24l) is valid 
only in a finite range of 02, restricted on both sides, and 
the width of this range turns out to be vanishingly small 
for ai ~ 0.5 - 0.8 used by MS. 

Specifically, Eq. (19) is obtained under three assump- 
tions: 

• that the T-dependence of the bosonic mass can 
be neglected and the bosonic propagator can be 



evaluated at T = - this is justifed if m sa 

• that typical vk are small compared to m, i..e., 
Dh{k) can be approximated by its small k form 
Db = v^k'^/{8TTm); 

• that typical k are small compared to fc^- 

We verified that the first assumption is well satisfied for 
all ai < 1 and arbitrary a2 and is not the subject of 
concern. 

To verify the second assumption we note that typical k 
leading to ([M)) are of order (mfci?J7„/w^)^/'^, and typical 



i iri 



T 



7=1/3 



for typical k 



Substituting the numbers, we find that 



vk 



2ain^/^ 



(25) 



where n > 1 is a number of a Matsubata frequency (typ- 
ical n = 0(1)). The condition z;A: << m is then satisfied 
when Qfi is small enough. There is no dependence on 
a2 in this equation, hence for sufficiently small ai the 
condition vk < m is satisfied for all a2- 

The third condition is satisfied at small a2 but breaks 
out at larger 012. Indeed, for typical fc. 



A ^ 27a2ni/3 



(26) 




FIG. 4: The interplay between Tp in the MS model and in 
the 7- model with the pairing kernel A(r2„) oc l/fij^. There is 
a very tiny range near 02 ~ 0.01 where Ti, in the MS model 



is the same as in 7 = 1/3 model: Tp = T^ 



IWEpai. 



At larger 02, the MS model becomes equivalent to the 7 
model with varying 7 > 1/3. The value of 7 increases 
with 02 and approaches 7 = 1 at large Q2. Tp in this 
regime deviates from ai dependence of Tp^ and eventu- 
ally saturates at Tp — T^^^ — 0.108-Ef. At smaller 02, the 
MS model becomes equivalent to the 7 model with varying 
7 < 1/3. Tp in this regime again deviates from qI depen- 
dence of Tp^ and tends to a non-zero value at a magnetic 



QCP, where 02 — 0. In this limit, 

equivalent to the 7-model with 7 = 

2 

n7=0+ ^. 1 o-Bp 



the MS model becomes 
= 0+ (A(n„) « logn„). 



and T„ = TJ 



1.. 



The condition k/kp << 1 is satisfied when the r.h.s. of 
is small. 



A. crossover at "large" 02 

The large numerical prefactor in (1^51) implies that Eq. 
(j24p becomes invalid beginning from quite small a2 ~ 
10^^. Once typical k become larger than kp, A(0„) has 
to be re-evaluated because the momentum integral in (fT9|) 
only extends up to kp. We assume and then verify that in 
this situation the Landau damping term J7„/A: is larger 
than Db(k) for all k < kp and for f2m comparable to 



the actual T„, which differs from T, 



7=1/3 



p 



in Eq. ([M)) . 



Neglecting Di,{k) in (|19p we obtain in this situation 



A(r!„) 



a 



7=1 



\nj 



(27) 



where fig 



4Ep/{3tt). Observe that J7q does not 

depend on 02, in distinction to fig" in Eq. (f^ . 

This form of A(ri„) corresponds to the 7— model with 
7 = 1. Tp in 7 = 1 model has been obtained in^ii^: 



T^=i = 0.254r22;"^ = O.lOS^i^. 



-.7=1/3 



scales as a^, the 



clearly see the difference: while Tp 
actual Tp saturates at large 02 . 

In the crossover regime, the coupling A(J7„) and the 
pairing instability temperature Tp should interplolate be- 
tween 7 = 1/3 and 7 = 1 behavior. The two tem- 



peratures Tp 



7=1/3 



and TJ"^ cross at 02 ^ 0.03, and 



the crossover should be around these a2- We show the 
crossover in Fig|31 Note that 012 ~ 0.03 is consistent 
with our earlier easimate of 012 ~ 10~^ at which typical 
k becomes comparable to kp. 

To understand the system behavior in the crossover 
regime, we computed A(r2„) for several 0:2 and found 
that for each given a2, A(J1„) can be well fitted over 
a wide frequency range by fl~^ form, where 7 changes 
between 1/3 and 1 when a2 increases. In other words, 
for a given a2, the MS model is equivalent to a 7— model 
with a particular exponent 7. In Figs . [5] and [51 we 
plot X{Qn) vs the number of the Matsubara frequency, 
and the actual Tp as a function of ai for a2 = 0.02 
and a2 — 0.06. The best fits to 7-model correspons to 



7 



0.6, 17: 



7=0.6 




0.04, and 7 = 0.7, fi, 



7=0.7 




0.11, 



respectively. The dashed lines in the figures for Tp are 
the values of the pairing instability temperature in the 



(28) 7-model (Ref.^^): 



0^7=0.6 
P 



o.5n 



7=0.6 




0.02£;f and 



This should be the actual Tp in the MS model at large 

a2. Comparing this TI''^"'^ with Tp ~ from Eq. (j24[) . we is some residual dependence of Tp on ai which is not 



^7=0.7 ^ 0.38172 = 0.042£:i.. We see that the agree- 
ment is nearly perfect, although for 012 = 0.02, there 
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FIG. 5: Solid lines: the dimensionless pairing kernel X{Qn) 
(panel (a)) for af = 0.32, and the critical temperature Tp 
vs the coupling ai (panel (b)) at 02 ~ 0.02 (02 measures the 
distance to an antiferromagnetic QCP). The dashed (red) line 
in panel (a) is X{Qm) = (f2o/n„)^. The best fit corresponds 
to 7 = 0.6 and fi^;""'^ Ri 0.045-Bf. The dotted (red) line 
in panel (b) is Tp in the 7 model for 7 = 0.6: T^=°-^ k: 
0.5r2g^ ■ = 0.02ii'F. The agreement is quite good, but there 
is some variation of the actual Tp/ Ep with a\ implying that 
the value of 7 in the fit of A(f2„) by 1/W slightly varies with 
ai. 



present in the 7— model. In Fig. [7]we show the fit to l/r2„ 
from of A„ and Tp for the limiting case of large 0:2 ■ We see 
that the actual Tp perfectly matches TJ^^ — 0.108i?F, 
and the ratio Tp/Ep is totally independent on ai. 



B. crossover at small 02 

We next consider what happens at small 0.2 , when typ- 
ical k are much smaller than kp- At the first glance, the 
analogy with 7 = 1/3 model should work in the small 02 
range because all three conditions used in the derivation 
of Eqs. ([2T|) - ((24l) are satisfied. If so, Tp should scale as 
a| and vanish at 012 = 0, see Eq. (|M|) . This, however, is 
not the case as is clearly seen from Figl3] ~ the actual Tp 
remains nonzero at a2 =0. 

The reason for this discrepancy is that the validity of 
the three conditions m » T,vk << to, and k << kp 
in fact only implies that Eq. ([M)) is self-consistent at 
small a2- However, there may be another contribution 
to Tp from the range vk > m, and if such additional 
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corresponds to 7 = 0.7, and f^Q 
(red) line is r;=°-^ ^ 0.38n2;=°-^ 



7=0.7 



FIG. 6: Same as in Fig. [5] but for a larger 0:2 ~ 0.06. The 
best fit by a power-law (dashed (red) line in panel (a)) now 

lAlEp. The dotted 
0M2Ef. The areement 
is quite good, and there is less variation of tp/ Ep than in Fig. 



contribution does exist, it may oversahadow the con- 
tribution from vk < m. A good example of such be- 
havior is the case of phonon superconductors at strong 
couphngiSiiiOe. McMillan T^*"^ - WDe-(i+^)/^ is self- 
consistently obtained as a contribution from fermions in 
a Fermi- liquid region of a; < ujp). Yet, there is another 
contribution to Tp from fermions outside of the Fermi- 
liquid regime, and at strong coupling (A >> 1) this sec- 
ond contribution overshadows the contribution from the 
Fermi liquid range and yields Allen-Dynes expressiorti^ 
T^d^ujdVX^T^'^'. 

To verify whether the same happens in our case, con- 
sider the limit 02 = 0, when Tp given by ()24p vanishes. 
If the actual Tp remains finite at this point, it neces- 
sary comes from vk > m simply because m{T = 0) = 0. 
Re-evaluating a(r2„) at ^2 = and T = we obtain 
Dbik) = vk/8 and 



A(f7„) = ^ 



l^«l.tanh- 



nvkp 




(29) 
At small fl„ << vkp this reduces, with logarithmical 
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FIG. 7: Dimensionless coupling constant \T{oJn) and critical 
temperature in the extreme limit of very large 02 = 600. 
Details are the same as the Fig. [5] The power-law exponent 
in the fit 7 = 1, f^^"^ = 4J5F/(37r) = QA2AEf, and rj7=^ f» 



0.254n; 

Qfl. 



,7=1 



O.lOSSF.This rj7= 



perfectly matches Tp for all 



accuracy, to 



K^n) 






^: 



7=0+ 




1^. 



where 51 



7=0+ 







(30) 



The model with such A(ri„) is again equivalent to the 
7- model, but with 7 = 0+. The 7 — 0+ model with 
logarithmical kernel has been analysed by Son^ and oth- 
It does have a pairing instability at a non-zero 



ersj 



.11.19 



temperature, i.e., the actual Tp does not vanish at a2 = 0, 
in distinction to Eq. (P^ . For ai — 0.6, we found 
TpT=°+ - O.OOSSj.. This obviously implies that Eq. ^ 
is valid only in limited range of a2 , bounded from both 
ends. Furthermore, compairing this T^^'^'^ ^ Q.005Ep 

with Tp~ from Eq. ([M)) . we find that they become 

comparable at the same 02 ~ 10~^, where Tp~ be- 
comes comparable to T^^^- As the consequence, the 
width of the range where Eq. ([M|) is valid turns out 
to be vanishingly small number-wise, although paramer- 
wise it is 0(1). This explains why the actual Tp is so 



different from T, 



7=1/3 



p 



see Fig. [3|). 



We will discuss the form of TIJ~'^~^ in more detail in 
the next section and here only note that the total Tp in 



the MS model is the sum of contributions from vk > m 

~ . In the crossover 



and vk < TO, i.e., Tp - T^=°+ 



■t;- 



range where Tj7^°^ and Tp~ are comparable, the MS 
model is most likely again fitted by the 7-model with 
varying 7 between and 1/3 (i.e., by different 7 for dif- 
ferent a2)- However, the width of the crossover range at 
small a2 is narrow, and we didn't attempt to fit the data 
at the smallest 012 by the 7— model with < 7 < 1/3. 



IV. THE MS MODEL AT THE SDW QCP 
EQUIVALENCE WITH COLOR 
SUPERCONDUCTIVITY 



THE 



The model with the logarithmical pairing kernel, as 
in Eq. ([50]) . is most known as a model for color super- 
conductivity of quarks due to the exchange by gluons. 
For small w„, the fermionic self-energy for logarithmical 
A(ri„) given by ((30)) is 



S(w„) 



duj'\{ijj') = 50; log 






(31) 



To shorten notations, we define g = ai/n . The pairing 
kernel at small frequencies is X{^m) — g log ( -^ — 1 . 

The equation for the pairing vertex is, to the logarith- 
mical accuracy. 



$(w) 



did' 



,Hoj')^Og{j^^)+\og{j^r^^) 



OJ' 



glog(H^) 



duj 



M-') ^°^^7k3 



l + fflog 



(32) 



The logarithm in the denominator is the contribution 
from the self-energy, Eq. ([31]). At small coupling, g « 
1, which we only consider, the self-energy can be dropped, 
and the condition that the solution of (15^ exists reduces 



more accurate solution was obtained by Son^: 



to g\og\^l-''^ /T^=''+) - 0(1), i.e., logT, ex l/v5- A 



TJ=o+ = c nr''+ e-^ =c^e' 
^ ai 



(33) 



where c and c = Sttc/c^ are numbers 0(1). The smaller is 
ai, the smaller is T^''^^"'" and the smaller is the crossover 

scale below which the actual Tp deviates from Tp~ . 



To obtain the prefactor c in ((33|) . one has to go be- 
yond the leading logarithmical accuracy and include (i) 
Y^ corrections to the argument of the exponent and (ii) 
the actual soft high-energy cutoff in the full expression 
for A(ri„), see Eq. ([29)) . The y^ corrections to the argu- 
ment of the exponent come from fermionic self-energy. 
These corrections have been computed analytically in 
Ref.^^, and the result is that the fermionic self-energy 
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FIG. 8: Solid lines: the pairing kernel X{^„) (panel (a)) 
and the critical temperature Tp in the MS model at the 
critical point, 0:2 = (panel (b)). In panel (a) we set 
Qfi = 0.32. The crossed (red) line in panel (a) is the fit to 



A(f7„) = (ai/7r^)log(f7;j-"+/n„) with ^ 



17 = 0+ 



SnEp 
— T 



(see 



Eq. H30|) V For comparison, the dotted (blue) line is the fit 
by a power-law with 7 — 0.6. Obviously, the fit by 7 = 0+ 
form is much better. The crossed (red) line in panel (b) cor- 

responds to T^=°+ f» l.Sf^e"^^^^. The fit is almost per- 
fect (solid and dashed lines are virtually undistinguishable). 
For comparison, we also show the fit to the BCS behavior 
Tp'^'^ tx l/aie""/^"! (the dashed (green) line), adjusting the 
prefactor to match Tp at ai = 0.75. We clearly see that Tp in 
the MS model follows T^=°+ rather than the BCS form. 



contributes the universal factor e^^^ +4)/8 — 0.177 to c. 
The effect of a soft high-energy cutoff depends on how 
the logarithmical form of A(ri„) is cut and is model- 
dependent. 

To the best of our knowledge, the non-trivial depen- 
dence of Tj/^^^ on the couphng g = ai/ir'^ has not been 
verified numerically. The equivalence between the MS 
model at q;2 = and the 7 model with 7 — 0+ allows 
us to do this and also to determine the value of the pref- 
actor in our case. In Fig. (|8]) we show the fit of A(ri„) 
for a2 = to l/r2„ dependence (which is almost perfect) 
and the numerical results for Tp/Ep vs ai. We fitted the 
data by Eq. ^, i.e. by (1/ai) e~''''^^'^^^'> dependence 
on ai, and, for comparison, by BCS dependence. We see 
that the fit by Eq. ([SS]) is nearly perfect. We emphasize 
that l/^yoi dependence in the exponent and the pres- 
ence of 1/ai in the overall factor are both relevant to the 
success of the fit. We view this agreement as the solid 



verification of Son's formula for color superconductivity. 
From the fit we also determined c « 0.037, and c « 1.8. 



V. CONCLUSIONS 

To conclude, in this paper we analysed the onset tem- 
perature Tp for the pairing in cuprate superconductors 
at small doping, when tendency towards antiferromag- 
netism is strong. We considered the model, introduced 
by Moon and Sachdev, which assumes that the Fermi 
surface retains pocket-like form, same as in an antifer- 
romagnetically ordered state, even when long-range an- 
tiferromagnetic order disappers. Within this model, the 
pairing between fermions is mediated by a gauge boson, 
whose propagator remains massless despite that niagnon 
dispersion acquires a finite gap in the magnetically disor- 
dered phase. From the point of view of quantum-critical 
phenomenon, the pairing problem then remains quantum 
critical in the sense that that the pairing is mediated by a 
gapless boson, and the same gapless boson that mediates 
pairing also accounts for the destruction of a Fermi-liquid 
behavior down to the smallest frequencies. 

We related the MS model to the generic 7— model of 
quantum-critical pairing with the pairing kernel A(il„) ex 
l/f2^ and the corresponding normal state self-energy 
S(a;„) ex w^"'''. We showed that, over some range of 
parameters, the pairing problem in the MS model in a 
paramegnetic phase predominantly comes from fermions 
with momenta k << kp, where kp is the Fermi momen- 
tum,, and vk << m, where v is the Fermi velocity, and 
m is the magnon gap. In this situation, the pairing in 
the MS model is equivalent to that in the 7-model with 
7 = 1/3 We found that Tp scales as a|, where 02 is 
the parameter which measures how far the system moves 
away from the SDW phase. 

On a more careful look, we found that the range of 
a2 where the analogy with the 7 — 1/3 model works is 
bounded at both small and "large" a2- At "large" 02 
(which numerically are still quite small), the condition 
k « kp breaks down. We demostrated that, in this 
situation, the MS model becomes equivalent to 7 model 
with varying 7, whose value increases with a2 and ap- 
proaches 7 = 1 at large a2- In the crossover regime, 
the actual Tp (0:2) deviates from the a| dependence and 
eventually saturates at large 0:2 at Tp ^ O.lEp. 

At the smallest a2, we found that the analogy with 
7 = 1/3 model again bveaks down, this time because 
the contribution to Tp from vk « m gets overshadowed 
by the contribution from vk » m. In this situation, 
the MS model becomes equivalent to the 7-model with 
7 < 1/3. Finally, at a2 — 0, i.e., at the SDW transition 
point, the MS model becomes equivalent to the 7 model 
with 7 = 0+, for which X(p.m) oc logfim- 

The 7 = 0-1- model has been studied in the context of 
color superconductivity and superconductivity near 3D 
ferro- and antiferromagnetic QCP. The pairing problem 
with the logarithmical pairing kernel is similar to the 



10 



BCS theory in the sense that the fermionic self-energy 
can be neglected in the leading logarithmic approxima- 
tion at weak coupling g. However, in distinction to the 
BCS theory, log Tc scales not as 1/g but as l/i/ff- We 
used the analogy between the MS and the 7 = 0-1- model 
at small frequencies and the fact that A(ri„) in the MS 
model does converge at high frequencies (i.e., one does 
not have to introduce a cutoff), and computed Tp nu- 
merically. We explicitly verified the l/y^ behavior in 
the exponent and also verified the 1/g dependence of the 
overall factor for Tp. To the best of our knowledge, this 
has not been done before. 

Overall, our results imply that the MS model is quite 



rich. It contains the physics of the quantum-critical 
pairing in all models with kernels A(f2„) oc 1/^2^ and 
< 7 < 1 . In addition, the MS model taken at the point 
of the magnetic instability is equivalent to the model 
for color superconductivity. It is amaizing that one can 
get useful information about color superconductivity by 
studying the pairing at the antiferromagnetic instability 
in the cuprates. We thank S. Sachdev, M. Metlitskii and 
J. Schmalian for the interest to this work and useful com- 
ments. This work was supported by NSF DMR-0757145 
and the Samsung Scholarship (E.G. Moon) and by NSF- 
DMR-0906953 (A. V. Ch.) 
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